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TECHNICAL NOTE 2869 


REFLECTION OF WEAK SHOCK WAVE FROM A 
BOUNDARY LAYER ALONG A FLAT PLATE. II - INTERACTION 
OF OBLIQUE SHOCK WAVE WITH A LAMINAR BOUNDARY LAYER 

ANALYZED BY DIFFERENTIAL-EQUATION METHOD 


By Yung-Huai Kuo 
SUMMARY 


By analogy with the boundary-layer concept, the flow produced by 
the interaction between a shock wave and a laminar boundary layer is 
subdivided into a viscous layer and a potential field. The assumptions 
that the compressibility effect in the inner layer is negligible and 
that the original flow in the outer layer is uniform lead to simple 
analytic solutions for the perturbed flow. The joining conditions at 
the interface between the layers determine an eigenvalue which gives 
the rate of decay and the character of the disturbances both upstream 
and downstream of the point of incidence. The final conclusions are 
in agreement with experiments. 


INTRODUCTION 


The present investigation is an independent study of the inter- 
action of an oblique shock with a laminar boundary layer in a compress- 
ible supersonic stream. In reference 1, where interaction of weak shock 
waves with both laminar and turbulent boundary layers was treated, the 
integrated momentum across the boundary layer was considered, rather 
than the balance among various dynamic effects at each point. This 
momentum-integral method is simple and, in certain respects, powerful 
and capable of yielding useful qualitative information such as the 
upstream pressure influence, pressure distribution, and the growth of 
boundary-layer thickness due to the presence of a shock, but it fails 
in regard to what actually happens inside the boundary layer. Іп the 
present report a different approach has been adopted, with the inten- 
tion of filling the gap left by the previous investigation. The pur- 
pose will, on the whole, be complementary, so as to provide a physical 
picture for the understanding of this complex phenomenon. | 
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Contrary to reference 1, the differential-equation method is 
employed here. According to available experimental observation, when 
an oblique shock is incident upon a laminar boundary layer the result- 
ant flow bears no resemblance to the flow predicted by potential theory. 
For if the viscous flow is absent the flow ahead of the shock will not 
be affected. Because of the presence of the boundary layer in which 
there is & subsonic layer, however, & sudden decrease of pressure at a 
point will immediately be transmitted forward by the inability of the 
subsonic layer to support an excess pressure rise. When the pressure 
is transmitted, the flow in the boundary layer will be retarded and the 
streamlines distorted. Since the outer field is supersonic, this change 
occurring in the viscous layer will affect the whole potential field. 
This is actually observed. For stronger shocks, the flow in the bound- 
ary layer generally will separate and will have backflow under the 
influence of an adverse pressure gradient. An adequate theory that is 
able to account for the observed effects cannot be formulated unless 
the boundary-layer hypothesis is abandoned entirely. One therefore is 
faced by a much more difficult mathematical problem. 


To restrict the scope and complexity of this study, let it be 
assumed first of all that the boundary is an insulated flat plate, and, 
secondly, that the incident shock is weak and its angle of incidence is 
such that regular reflection would be possible, were the flow friction- 
less, and, lastly, that the free-stream Mach number is not large. Under 
these assumptions, the angle of deflection of the flow in passing through 
the shock wave will be small, and the temperature variation between the 
free-stream condition and that of the plate will not be large. In fact, 
the study of a laminar boundary layer indicates that, for moderate Mach 
number, the temperature as well as the compressibility effects are unim- 
portant (reference 2). This must remain true even if the flow is not 
boundary-layer flow. Therefore, without loss of generality, it will be 
assumed that the viscosity and thermal conductivity of the gas will be 
taken as constant and the Prandtl number is unity. 


In order to bring the interaction problem within the scope of prac- 
tical mathematical analysis, these simplifying hypotheses have to be 
made in the absence of a proper method of approximation, such as the 
boundary-layer theory. Broadly speaking, examination of schlieren 
photographs of the flow produced by the interaction between a shock and 
a laminar boundary layer will reveal that two characteristically dif- 
ferent outer and inner regions exist for sufficiently high Reynolds 
number. Тһе outer field is characterized by its strong potential char- 
acter, whereas the region close to the wall is predominantly viscous, 
which is quite reminiscent of the boundary-layer flow. It appears 
natural, therefore, to assume & priori that the viscous effect is con- 
fined to a thin layer in the vicinity of the boundary and the outer main 
flow is potential. These two different flows are then in dynamical 
equilibrium. If one is disturbed, the other will be affected. Since 
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the outer field is supersonic, any local change will be felt in a much 
larger region then that in & subsonic field. 


After the flow field is separated into two regions, specific assump- 
tions regarding the structure of the viscous layer can be made. It is 
important to note that in the case of a compression incident wave, the 
overwhelming effect taking place in the viscous layer is the sudden 
decrease of the velocity or even reversal of flow (backflow). If back- 
flow sets in, the flow speed in the subsonic region will be very much 
reduced. As a result of this, the streamlines will be pushed outward 
and the flow compressed. On account of tbe displacement of streamlines, 
the subsonic region will become thicker, and the thickening of the sub- 
sonic region is & decisive factor that distinguishes the strong from 
weak interactions. The importance of this dimension (the thickness of 
the subsonic layer) has already been established by Tsien and Finston 
(reference 3) in the case of inviscid theory. 


The viscous layer, thus, must have two distinct sublayers, each 
displaying a different character. In the supersonic layer, the flow is 
characterized by large velocities and is not unlike that of an unsepa- 
rated boundary layer. Therefore, for this layer both viscous and inertia 
forces.are important. On the other hand, in the subsonic layer, espe- 
cially with backflow, the average speed will be very small. Tn this 
case, because of slow motion and moderate temperature change, the change 
Of density is always a lower-order effect. In fact, the contribution 
due to compressibility is proportional to the square of Mach number, 
and, if the average Mach number in the subsonic region is small, the 
compressibility effect is, indeed, negligible. Because of this physical 
fact, the subsonic layer will be taken as incompressible. 


With these assumptions, the problem is finally solved by pertur- 
bation of weak incident waves. As a test of these assumptions, a simple 
flow with broken-line velocity profile is taken as the basic flow: In 
the incompressible layer, the velocity is a linear function of the dis- 
tance from the plate; in the compressible layer, it is constant. The 
density in the basic flow is constant in each layer but takes different 
values. Thus, at the interface where the two flows join, the velocity 
is continuous but density is discontinuous. For this case, a first-order 
solution consistent with these assumptions is completely determined. 


In the case of weak shock, excellent agreement with experiments has 
been achieved for the pressure distribution on the plate. It confirms 
the conjecture that separation of the flow as well as backflow always 
occur. Because of the occurrence of backflow, transition downstream of 
the point of incidence is unavoidable in the viscous layer. There are 
strong experimental evidences but detailed investigations are yet to 
be conducted. In the outer field, on the other hand, it is predicted 
that in the place of the regularly reflected shock there is a strong 
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expansion, and farther downstream & train of strong compression waves 
must exist, eventually forming an.envelope. Therefore, downstream of 
the point of incidence а second shock must occur. This deduction is 
also confirmed by experiments. 


Lastly, the importance of nonlinear effect is discussed. 


This study was conducted at the Graduate School of Aeronautical 
Engineering, Cornell University, under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
The author wishes to acknowledge the efficient assistance of Mr. Nelson 


H. Кешр. 


SYMBOLS 
А constant 
а! speed of sound 
Сч, C2, 05, C, C5, 06 constants 
D constant 
d constant 
а/аъ convective derivative 
F,G scalar functions 
figi defined by equations (43) and (41), respectively 


H nondimensional enthalpy 
; 3/8 
11/3 ЕШ Bessel function of first kind with imaginary argument 


J1/3 Bessel function of first kind with real argument 


KJ Ko , Ko, K), constants 


NACA TN 2869 S 


3/2 
TE ) Bessel function of second kind with imaginary argument 


Ge г i : Pe 


Қ = 
L length of plate. from leading edge to point of incidence 
M Mach number 
n dilatation 
P pressure 
Pi, Poy Р, constants 
p nondimensional pressure (s Бл) 
p' pressure 
1⁄2 
R Reynolds number (U L/v,) / 
2 
rx r(2/3)/3 [3 
T nondimensional temperature 
T temperature 
8,t defined by equations (Bh) 
U velocity 
|u,V nondimensional velocity components 
u',v' veloeity components 
(1) 1 
U `, v constants 


Б dne = - — ae M — ilt colin i nn — os 


a 
! 
= 

Ww 

W 
Do 


аг 
| 
o 


> 
il 


nondimensional Cartesian coordinates 
coordinates 


defined by equation (Ch) 


üefined by equation (30) | 


defined by equation (83) 


defined by equation (C15) 


gamma function 

ratio of specific heats 
Laplacian operator 

deflection of flow (equation (B6)) 
flow-deflection angle 

defined by equation (31) 
defined by equation (C6) 
defined by equation (C5) 
eigenvalue 

free-stream kinematic viscosity 
nondimensional density 


density 
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= XR - XM 


1/3 


а 
І 


+4 
l 


= AR) 


(D 
Subscripts: 
C 


с 


density at у = 0 


defined by equation (CT) 
nondimensional velocity potential 


stream function 


angle between velocity and shock 


complementary 

solution in compressible layer 
solution in incompressible layer 
value on plate 

particular 

due to step wave 

at separation point 

due to transmitted wave 

partial derivative with respect to 


partial derivative with respect фо 


x 


J 


in region 1, before shock (see fig. 1) 


in region 4, behind shock (see fig. 1) 


in region 3 (see fig. 1) ) 


free stream I 


———— — — ре — — a = ee ш — —— ------- 
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|| parallel to oblique shock wave " 
ab perpendicular to oblique shock wave 

Superscripts: 

(0) zeroth order 

(1) first order 

(2) second order 

(3) third order 


STATEMENT OF PROBLEM AND BASIC ASSUMPTIONS 


Let there be & laminar boundary layer in & compressible viscous 
fluid along an insulated flat plate immersed in a steady uniform super- 
sonic stream, and let an oblique shock be incident upon the plate. When 
the steady condition is established, the flow in the neighborhood exhibits 
а character which is entirely different from the original flow. 


From all indications, this flow does not obey the boundary-layer * 
approximation; nevertheless, for simplicity for future analysis, the 
concept of boundary layer, or viscous layer, will be retained.  Namely, 
the whole flow field is visualized as consisting of inviscid and viscous 
flow in equilibrium with each other. The possibility of existence of 
such a demarcation line will be assumed at this moment. The solutions 
obtained are consistent with the assumption, as will be seen later, so 
that the theory is self-consistent at least. Naturally, its further 
justification rests upon experimental evidence. 


The main feature in the viscous layer is that, in the case of a 
compression wave in the inviscid outer flow, backflow generally exists. 
For this reason, terms in the equations of motion which are unimportant 
according to boundary-layer approximation become decisive as the sup- 
posed large-order terms vanish. Therefore, the pressure gradients along 
both directions have to be considered. То simplify the mathematical 
process, some minor effects, such as the variation of the viscosity 
coefficient and thermal conductivity with temperature, will be neg- 
lected and the Prandtl number will be taken equal to 1. For moderate 
Mach numbers, less than 3, say, this neglect, according to boundary- 
layer investigations, will have little effect on the major character- 
istics.of the flow. 


The inviscid flow generally is rotational, as it involves shocks. 
This is particularly true in the case of & local supersonic zone over 


AE 
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a curved surface. For a purely supersonic flow, if the shock is 
slightly curved, the vorticity generated behind the shock is of high 
order and can be neglected. The perturbed flow in the outer region can 
then be regarded &s irrotational. 


Ав the region which is influenced by the presence of the shock, 
according to experimental observations, is confined to айп area about 
the point of incidence, with a dimension only a fraction of the length 
of the plate, a point which lies at в distance about the length of the 
plate from the point of incidence will be considered as at infinity. 
Consequently, the boundary-layer flow will be replaced by a "shear flow" 
extending to both positive and negative infinities. This approximation 
is justified if the derivative along the flow is much larger in the 
perturbed flow than that in the original flow. For large Reynolds num- 
ber, this condition can always be satisfied. 


METHOD OF SOLUTTON 


The flow is supposed to be two dimensional and steady and the fluid 
is compressible. If the flow field can be subdivided into viscous and 
nonviscous regions, the flow in the viscous layer satisfies in dimension- 
less variables the system: 


р ші = -Px + gn + 5 nx) (1) 
оз” Py + (Дт +з ny) (2) 
(ou), + (ov), = 0 (3) 
ан 1 (y - 1)мё |/au àv 2 2 
o s gar + SH) t к= уу үн 


and іп the case of perfect gas 


p = (ae) m (5) 


_ - 3 -- —~ -- aw — p zm жез. oe 
-- ------  - - -— —— —— — — - Vn- Ó—— —— — ——[ n: —À —R—— — a 
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Here the subscripts denote partial differentiations with respect to the s 
Cartesian coordinate x or y indicated; the Laplacian operator A, 
the convective derivative d/dt, the nondimensional enthalpy Н, and 
the dilatation n are defined by 


39 
бу (6) 
м 


Furthermore, the velocity components u' and v', the pressure р', 
the density р!, the temperature Т', and the coordinates x' and y' 
are nondimensionalized as follows: 


g = n" 
U, : 
m ли 
Ue 
1 
р--Е- 
MM 
t 
p = E (T) 
Poo 
х =-® 
Vol 
Us 
t 
y = — 
Vool 
Оо 
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where Us, р, and v, are free-stream velocity, density, and kinematic 


viscosity, and L is the length of the plate from the leading edge to / 
1/2 

01 

the point of incidence. Finally, the parameters y, М, and R= TE 


со 


stand for, respectively, the ratio of specific heats, the free-stream 
Mach number, and the Reynolds number. 
On the other hand, tbe inviscid flow in the outer region is assumed 


to be irrotational amd, consequently, is determined by a nondimensional 
velocity potential Ф(х,у) satisfying the equation 


[Ж : Y les - ачуу + (еу à Y les, = 0 (8) 


where the sound speed а! is given in terms of the velocity components 
п and у! by the relation 


(a5 = (au!) - — ty + (1) - ve] (9) 


By the assumption that these two different flows are in equilibrium with 
each other, the flow determined from equation (9) must join smoothly with 


that given by the system of equations (1) to (5) subject to the boundary 
conditions | 


when y = О and 


uz 
у = Фу 


when jy = o. 

As the disturbances are initiated by the incident shock, it is 
expected that they are small if the incident wave is weak. From experi- 
ence, this is &t least the case in the field upstream of the point of 
incidence. Т? the shock strength is characterized by the flow-deflection 
angle є, then for small values of є the solutions are expanded in 
powers of в: 

Pd 
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u = ul) + eall) + eal) 6... 


v 2 v(0) + evll) + 42,0) 4... 
p = р(0) + ep (1) 4 cep (2) +... (10) 
р = o (9) t eo) + є2р(2) +... 


Т = т(0) + ет(2) + eO! к... 


By substituting the quantities defined in equations (10) іпбо the system 
of equations (1) to (5), there results the following system of equations, 
according to the powers ОҒ s; namely, for the zero order, 


(0) (0) | 
0) мы 0), Қы + x 


5 (0) E i р, (0) " ils) + ła (0) 


3 y 


(ь(ом(о)) 4 (s(0,(0)) = o 
| | (11) 


dt dt 
x y 


о(0) а(00(0) -1лн(0), Z- 1)м (Еш x (49,9) _ 
dt R R 


2(.(0)4(0)) ем 
3 х 3 y 


ь(0) = (au) L5 (Og (0) 


— a a — —- —— — —— — — — F —— x 2. —— — — a 
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where 


2 2 (12) 
10). (0), 223 2.07, (oy) 
Csa (9), (0) 
for the first order, 
Oa AC? А Oo) LA EUNT б}, 
d dt dt a 
Қай, 2 3,07) 
purs aD DON, sk, 
d dt dt J 
Loe) ‚1 a 
(13) 


(00202) + pul) + (2000) + 2000), = o 


CTO K (ara , 0) Jam „з шш), 


(y - 1)MÉ (49,0 E КЕ EE 
R dt ау À | а à 5 


2(y(0)n(2) + w(t)y(0)) | _ 2 (o) (1) n exe) Я 


02) = (эё) (KOK) + ,G)x(o)) 


14 


where 


a) уз, ә 
dt ze 


wD) 2-0), 6 2) ОМ). (бй) 
a(t) 2,2) 4 y (1) 


and for the second order, 


(0,00 4,QXQ) ,(0)Q) (0), (1) 
(Est, em OP), ы 


dt 


А Sue uid 00) a(0), (0) 
d dt dt 


-p 0) + (62) " 1 n, (2)) 


oe аба), e - aae e? | «e M (0,0) _ 


dt dt dt 


FROM Har?) : 2 ny (2 ) 
(200062) 4 kG), р(20200)) + (2000,62) 4 (13,0) + (2 (о), m 


0) 


(0), (2) (2) ( (2) (1) (0) 
p (9) Cue Sid f К? a + p(t) 40), G) n К a + p(t) ар S) 0). to) = 


dt 


paal) , (z - ay (49%) „030, 60) . 


dt x 


(0) (2) (1) (1) (2) ud 
dV ^v d^. d^ ^. 2( (0) (2) (1) (1) (2) 
( dt + — ^ + LE ен * - 2, n + u n Tu а) - 


2(,(0),(2) + y 0,0) , yl КОҢ 


pl) „ (ng) ^ (, (05) , ,Q50) , ,(2)x(o)) 
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(14) 


(15) 
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where 


a2) суа, (2) ә 
4% Ox ду 


2 2 
н(2) = т(2) 4 (у - TN 3:940). 2) ; ай, 07. (16) 


202) (8) 4 y (2) 
In the case of the potential flow, if the velocity potential is 
expanded as | 
Q = x + өр\ї) + gag (2) к... (17) 


there will occur, similarly, equations for the first- and second-order 
quantities. These are: For the first order 


2 (1) (1) zd (18) 


and for the second order 


NONO EE OM s 0,0] (19) 


where В = Y - 1. 


FIRST-ORDER SOLUTION - UPSTREAM, x <0 


Iet the point of incidence be chosen as the origin of the Cartesian 
coordinates. Then, negative x will be called upstream, and positive x, 
downstream, of the point of incidence. The various regions will be num- 
bered by 1, 2, and 3 as shown in figure 1 in which OS indicates shock 
and OM, the limiting Mach line of region 2. 
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Ав it has been assumed that the basic flow is a laminar boundary- 
layer flow of a Blasius type, such as considered by Von Kármán and Tsien 
(reference 2) and by Emmons and Brainerd (references | ала 5), system (11) 
simply reduces to the well-known Prandtl boundary-layer equations, whose. 
Solution has already been found in these references. 


Since the basic flow is & boundary-layer flow with constant pres- 


sure »(9), system (13) can be simplified. From experiments it has been 
established that tbe space variation of the perturbed quantities is much 
more rapid than that of the unperturbed quantities. By the boundary- 
layer approximation, the ratio of the partial derivatives 0/dx:d/dy 

is 1:R. For large Reynolds number, the x-derivatives of the basic 
flow, as well as у(0) which is of the order В-і, can be neglected in 
comparison with the y-derivatives. That is, u (0) and p 9) are func- 
tions of y only. This approximation is confirmed by experiments and, 
as а matter of fact, it is customarily used in the pressure measurements, 
because the wave is shifted forward and back relative to & fixed pressure 
orifice to measure the pressure distribution before and behind the wave. 


In the case of unit Prandtl number, g (0) 
dropping terms such as y (044,01) and ux (O) t), system‘ (13) becomes 


is a constant. Then, by 


(ula) + u (0)()) = >, G) + Ma) + 2 a, G) 
(009), (0) = py) + Howl) + L ny(2)) 


(0) + wp (1) + 02) (0) = o 


(20) 
ҚОҚ) (2) = 1 лн) + aM o) а) + vy O 
(9, (1)) - Zan) | 


У 


ph) = (ag) 000,090) 
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Incompressible Layer 


It is noted that the coefficients of system (20) depend on the 


basic velocity profile s Since (0) (у) cannot be expressed 


by simple functions, in order to simplify the analytical work further 
simplifications are necessary. 


First of all, in the case of insulated plate, the temperature and, 
hence, the density in the basic flow will have a vanishing gradient on 
the plate. This makes the variation of the density in the viscous layer 
much smaller than that of the velocity. When the flow is subject to an 
adverse pressure gradient, the flow will be further retarded. When the 
backflow occurs, the average speed will be very low and, consequently, 
the representative average Mach number will be small. Under this condi- 
tion, the change of density is less important than that of the pressure. 
A layer where this approximation is valid is called an "incompressible 
layer" and the dominant effects will be pressure and frictional forces. 


Of course, it is difficult to define the thickness of this layer 
beforehand. Generally, it would correspond to the subsonic portion of 
the viscous layer. Іп the basic flow, the sonic boundary сап be exactly 
calculated. When the flow is perturbed, it is unknown, but certainly 
will be thicker, for the flow is subject to an adverse pressure gra- 
dient. Owing to this fact and also for mathematical expediency, the 
thickness will be taken in accordance with the way the velocity profile 


u(O)(y) is represented. "This will become clear below and the thickness 
defined turns out to be nearly half the original boundary-layer thickness. 


Because of the assumption that density is constant and assumes the 
value o, (0) on the plate, system (20) simplifies to 


o, C) (ula, (2) į u (O)y(2)) = -p U) + &ли(1) 
pul, (2) = ap (2) + 1 ala) 
(21) 


n QJ) = 0 


5040 (1) = а a (ü) , 20 - Н (о), (1) 
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Thus, the velocities, pressure, and temperature can be dealt with inde- 
pendently. For the velocities, tbe first three equations yield by 


elimination of р); 


o OO - 2,00) о (0%, (040) = Dafa) - + G)) 
(22) 
5 Q) = О 


After the velocities are known, the temperature will be given by the 
energy equation. Іп the present problem, however, temperature is not 
am interesting quantity and will subsequently not be mentioned. То 


solve equations (22), let а stream function (3 y) be introduced 
by the relations 


a(t) y 0) 


Eo a (5 
From equations (22) y then satisfies the equation 
о, (90.00), (2) = BA Ау) + p (Ou, (0, 1) TT 


By neglecting the compressibility effect in the inner layer, the 


basic velocity profile u(9) (y) becomes the Blasius profile. The vari- 
able coefficients of equation (23), for small values of y, will be 
power series in y. However, it has been recognized that in the case 
of Blasius profile the velocity attains the free-stream value fairly 
rapidly and the initial portion is nearly linear. Consequently, it will 
not involve serious error to replace the continuous profile by & broken- 
line one, so that ihe initial part is proportional to y and the 
remaining part, constant with the free-stream value unity. If the skin 
friction agrees with the exact value, the velocity profile can then be 
defined as 


02705 = су when OŠ у < Уу 
(2%) 


ul) (у) =1 when y, < y < o 


”“ 


a m “-.---- 2 = — = | — -- — —— — — — —PyaR-x a —F -- --- - — ———— - ——— — — 
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If the velocity is continuous at y - Ji» then yi = егі which is about 
half the boundary-layer thickness, namely y, ^ 3. TF the incompressible 


layer is defined as the interval 0 Sy 5 y, then equation (23) becomes 


бу Ay, V) = д ay? | (25) 


where dà = par. 


To solve equation (25), it is natural to assume the form 


ytt) = ү(у)е^Х 
(26) 
һ > О 


for region l where x < 0. This assumes a special form of compression 
waves induced in the potential field; that is, the solution of equa- 
tion (18) 


p(t) E А «“х-Ву) . (21) 


where А is a constant to be determined. This is entirely in agreement 
with the conclusions reached in reference 1. Substituting equations (26) 
in equation (25) and simultaneously putting 


Vo +r = z(y) (28) 
there results 
z + (лау - 2), = О (29) 


The general solution of this equation is 


- TAE 13/2) И сал/ 2/2 15/2) (30) 
3 3 


— — — —  ——.— cay Rn — n  — —  — m шел -— — — —— 
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where 
7 = Ga.) Xy - Xx 1) (31) 


C. and C, are constants, and 12 43/ 4 and x, (2 id^ е denote the 


Bessel functions of the first and second kind: with imaginary argument. 
With z known, the solution of equation (28) is 
i 


--— 5 


V = C3 cos Th + C, sin тд + 


in 7 | T 

SETA [ (eos 0G an - SETI f (sin тп) (а) ап (20) 
No No 

with т = x63) 1/3, On the plate y = O and 402) = y Q) = 0; hence 

Ұ(0) = ү'(0) = О. For large Reynolds number and small values of А, 


y = О corresponds approximately to тп = 0. Тһе boundary conditions 
thus require C3 = Cj = O. Moreover, at у = Ji» 1 will be large, 


7 COS 
and since тту = Ay,« 1, the integrals J 1/2 TT) ii (2 73/2) dn 
O sin 3 3 


2 „3/2 
E = 
will.diverge like J e3 1 ап ав т approaches infinity. Therefore, 
0 


Co = 0. The solution in the incompressible layer can then be read ав 
follows 


T 


елес TUE EET Иа PEN 
y (s п) | к (ес oe (d 9/2) ал 


n 
cos т, 1/2 (sin тт ЗЕ 13/2) ап [е^ (32) 
; 3 
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Thereby the velocity components аге 


404) zi с. (вә) /3 к= UA 1/2 (cos "2 15/2) ат +. 
3 


т 
(sin Tq) I 1^ (stu т) (2 73/2) ne 
O. 3 


t (33) 


“- 


(H) -= -c, (ду.)1/3 Z 5) IÑ (cos "2 13/2) dq - 
3 


(cos тт) ГА 11/2 (sin Tm x 73/2) чрн 
3 


On the other hand, the pressure can be obtained by integration of 
the first of equations (21) and is, by the form of solution, 


ou, (02,0) . 
X 


p(t) = 9 (0 (0,0) _ 2 СЫ CURT 


where the constant of integration vanishes by the condition &t negative 


infinity. Substituting 402) and V) from equ&tion (33) апа making 


use of equation (28), & straightforward reduction yields 


(Ос n 
1 1/2 2 „3/2 
COS T (sin ттү)К-| & ад - 
(con m) / mG «ж 6 03/2) as 


TT (ее f^ 72/2 (сов т) [Š 73/2) йт + 
3 


(sin Tm) J. «vM (stu та) [2 93/2) 3 ex (3%) 
3 


—— - —— ———ьһ 
— —— - = e ——— e — —— Y. r nA... 


— - —  — — = — с —— P P P  — - 
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On the plate y = n = O this reduces to 


» 06,0) = =, acy ie ($ 3/2) ех 
3 
=0 


Compressible layer 
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(За) 


In the compressible layer, because of the forward momentum of the 
potential flow outside, the velocity is everywhere positive and differs 
from the basic flow by & small amount in order to support the pressure 
rise. As the velocity is high, both compressibility and viscous effects 
will become important. Under this condition, the problem is considerably 


simplified by taking a uniform basic flow, namely 60) = 
discussed above. Accordingly, equations (20) reduce to 


а, V x -p V) T гч) + 3 a.) 


vy = (2) + (1) +2 ny) 


ко 45 Qe 


тр, 


(1)-. 


н„ (1) 1 an) , & zu» MO 


pilos Ca В Қо) 


where 


нї) at (y - iy RO? 


(0) 


= 1, ав 


(35) 


The elimination of р(1) from the first two equations gives 


[Ú G) _ v GJ) - 1 afu) _ S) 


x 


(36) 


=- —— — — ——— ------ 
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APter differentiating the first and the second of equations (35) with 
respect to x and y, respectively, an addition yields 


0). („4 мб) 


moreover, by definition 
g, VJ Е yp, (2) + (у - 1)ма, C) pnt?) 


A substitution of a(t) апа ли, (2) in the energy equation leads to 


- wu (1) + э COME a) - Lan) қ а Mn, (1) + 


и LEPP ANG) зт) 


ху 


Equations (36) and (37) form a system of linear partial differential 


equations for y(t) and vit with constant coefficients. The solu- 
tions, if carried out, are expressible in terms of trigonometric and 
exponential functions. According to the arguments of ‘these functions, 
they fall into two groups: One varies slowly and is "potential-like"; 
while the other varies rapidly and, therefore, is "boundary-layer-like." 
The latter group consists of two exponentials whose arguments differ by 


& term of о(2). Hence, for large Reynolds number, the latter two expo- 


nentials will degenerate into one. It was found that this form of solu- 
tion can be obtained by solving а much simpler problem, namely, by 


assuming y(t) = 0. This assumption appears to be nothing but a method 


of approximating the solution of equations (36) and (37) in the case of 
large Reynolds number. 


re g) 
is in its place 


is taken to be zero, then instead of equation (37) there 


ES + мъ (1) = LN $3 ЖЫ) (38) 
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by elimination of pt) and ml), To solve the system (36) and (38), 


two scalar functions p(t) and g(t) are introduced through the 
relations 


ult) = p (1) 4 ay (1) 


(39) 
т lg @) 5,09 
The equivalent system, then, is as follows: 
1 1 
ag, V = lA ag! 
(40 ) 


According to the form of the potential solution, it will again be 
assumed that 


By the condition that (2) is finite at infinity, the solution is shown 
to be 


g = C3 сов Xy + Су sin Ay + Be (41) 


where Со, C3» and С, are constants of integration and 


of = AR - А (42 ) 
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By substituting Бу from equation (41) in the second of equations (NO) 
the solution for fi was found to be 


fy = C3 sin Ay - Ch COS Xy + ey + Cee 


(у - 1)м 90 oy 


ee (43 ) 
À 
х + (F з) 
where 
TN. | yM, а. и 
l \ 3 3R 
The velocity components for у> уу, by equation (39), are given by 
woe 3) «acf PUn) " 
а (а - >= 9 sce 79-71) (hh) 


s e 
vil) 2 


— L — - 1)м Мо. 0 71) ex 
2 
XR «( (52 : 1 MR 


- pis Mr) + аде 15971) қ 
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by redefining the constants С. It is noticed that the first group of 
exponentials varies with both x and y with & slope proportional 

to А, which, according to experimental evidence for this type of flow, 
is a small number; that is, the variation with x is relatively slow. 
On the other hand, the second exponential varies with y with deriva- 


tive proportional to о, which is of o( R). For large Reynolds num- 
ber, it appears that 0 could be thousands of times larger tham 2. 
Therefore, in the case of large Reynolds number where the viscous layer 
is thin, the first group of terms changes only slightly while the second 
exponential drops to zero. Applying the boundary-layer concept, the 
potential-like terms will be taken as constant and equal to the boundary 
value of the potential solution. Consequently, the solution can be 
written as 


(Ма) 


with 


It is seen that when o(y - Yi) becomes large, this solution joins the 


potential solution at' y = O and the constants Cs and Сс мі11 Те 
considered as eliminated. Ву means of this approximate solution the 


pressure р(1) is shown to be 


Me» = fs - *X ше 678 ue (45 ) 
where 
ae x - 2 зае) 
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Eigenvalues X 


In the last two sections, the velocity and pressure have been cal- 
culated by entirely different methods of approximation in two different 
layers. For the incompressible layer, the nonslip conditions are sat- 
isfied on the plate, whereas in the compressible layer, the velocity 
agrees with the potential flow at infinity. The complete solution is 
then left with three undetermined constants Cj, Со, and A and with 


an arbitrary parameter А. To determine these constants, it is assumed 
that the two solutions (33) and (34) and (hh) and (45) must join at the 
interface y = Уу: Now, because of the simplification made in connection 


with equation (lh), the conditions at the interface are to be stated as 
follows for y = yi: 


Қорғас 
v(t) = v Q0) (46) 
92,9 


where the subscripts с and i indicate, respectively, solutions of 
compressible and incompressible layers. Тһе u velocity profile thus 
will have a discontinuity in slope. This could be improved by dropping 


the assumption н(2) = 0. However, &s pressure depends only on the 


velocity and its second derivatives, an error in shear can produce only 
minor contributions and can be ignored. Ву substituting the solutions 
in equations (46), there results the following system of equations: 


-BA - (xs? + т\с + ves yo = 0 (47) 


~, Ко (1) Е 
-KA + — Ca = P (п) = 0 
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where U), үй) ana pl) аге definea by 


(1) (ga 1/3 сов T п 1⁄2 s T))K, d 
U (xr) (cos w)J. 1 (cos TH 1 п + 


11. 
(sin m), "ШУ? (sin TN Жі dn 
3 


yt) Е (а) /3 (віп т) ^ 21/2 (cos TH Жү ат - 
3 


7 
(cos ту) f. d 11/2 (sin тп Жү ат 
3 


XO 


T 


PG) 


MEE n 3/2) РА (ay 72/3 y (2) _ тт, (ах) 300) 


In order that linear system (47) will admit a nontrivial solution for 
Сі» Со, and A, it is both necessary and sufficient that the determinant 


Deu om 


B r.a) ү) = 0 (48) 


— - ——— —— — YY — -- - 
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vanish. This equation serves to determine all proper values of А. 
After AX has been determined, the constants Су and Со corresponding 


to this А can be solved for in terms of А, namely 


C4 = A ——— T 
и - (aE + go 
(49) 
| Е. putt) _ y(2) 
(e ci na 


where the constent А remains to be determined. 


To solve determinantal equation (48) analytically even for this 
simplest case does not seem to be possible. The procedure from here on 
is essentially numerical. For the present purpose, the numerical solu- 
tion will be carried out, based on the approximation that 74 will be 


taken mathematically as infinity and TT, small. Accordingly, the 


integrals in y(t) v QJ, and р(1) will have ап infinite upper limit 
and cos Th, 5 1 and sin тту % ХУ). Furthermore, for large values 


" -3/% -2 113/2 
n e e 


of 14 the Bessel function «(2 will be 


3 


very small and can, therefore, be neglected. Equation (h8) 1s finally 
expressed explicitly in terms of the parameters R, p, and a. То 


retain terms up to the order olr 8/6), the determinantal equation 
becomes 


where г- г(2)/32/3. 


For M= 2, R= 774, and a = 0.332, there are two pairs of posi- 
tive and negative roots. The negative roots would make the perturbation 
infinite at negative infinity, which is contrary to assumption. Therefore, 
negative roots are not admitted. For the pair of positive roots, one 
is roughly 20 times the other which is 0.0467. For large Reynolds num- 


. .-3/8 31/2 зі -1/h 
ber, these would correspond to. X4 В 3/ R / / R / š 


ry a3 + À + Жы: = 0 (50) 


and л B 


— — — —M _— “m — 
ee Cm — 
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In region 1, x is negative. The disturbance with a logarithmic rate 
of decay A, will quickly disappear. The observed disturbance must 


have a decay rate equal to Хә. 


According to this solution, the dependence of the distance influ- 
enced by pressure disturbance on the Mach number and Reynolds number can 
be discussed. By the form of the first-order solution, the pressure 
decays exponentially from the point of incidence. For & given value 
ОҒ е, when -АХ = d, d being a constant, the pressure disturbance 


would have dropped to а certain fraction of its initial value. ‘Thus, 
-X = а/Хә would serve ав а measure of the distance reached by the pres- 


sure. Therefore, by varying M and R, the distance X'/L will change 


according to the law p 3/g-3/^. Namely, by increasing both Mach num- 
ber and Reynolds number, the distance reached by pressure disturbances 
will decrease. This result is quantitatively different from that arrived 


at in reference 1, which is g 28-12. By comparison there is an 
increase of both compressibility and viscous effects in the new result. 
This shows how gross an error can be made if the boundary-layer арргохі- 
mation is applied. Tt is surprising also that the upstream pressure 
propagation depends only on M and R but, to the first order at least, 
is independent of the shock strength. This seems to be in agreement with 
experiment (reference 6). 


It сап be shown from equetions (33) and (34) that for arbitrary 
values of А the pressure р (10) at the point of separation 


where цу = 0 is proportional to в71/2р-1/2, 
FIRST-ORDER SOLUTION - DOWNSTREAM x > O 


As in the upstream first-order solution, one must start with the 
potential solution. In order to determine the flow in region 2, the 


interaction between the incoming wave o! and the shock must be 


considered. 
Potential Solution 


In equation (27) the incoming wave is given by 


190) = ^ JA (x-By) 
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When this train of waves hits the shock, the shock will be slightly 
modified according to Rankine-Hugoniot conditions. The velocity poten- 
tial in region 2 with these conditions satisfied is shown to be (appen- 
dixes A and B) 


(1) 


QU cg + By) + 2 eh (x-by) (51) 


X 


That is, to the first order, the flow in region 2 is simply & super- 
position of & step wave upon the transmitted wave. Since the conditions 
are prescribed on the shock 05, they are uniquely defined in region 2, 
terminated only by the extreme Mach line OM. 


To continue this solution into region 3, two possibilities present 
themselves. It шау be either continuous or discontinuous there. If it 
were continuous on ОМ, there would be а discontinuity in pressure аб 
the point of incidence. Іп the case of inviscid flow, this, of course, 
would be admissible. But, by the condition that the viscous layer will 
Join smoothly with the potential field, this would make the velocity 


uc у) discontinuous at the origin x = 0. Hence, the first possi- 
bility must be discarded. If the solution is discontinuous at ОМ, and 


the pressure as well as the velocity "n &re continuous &t the origin, 
then the discontinuity must correspond to an expansion. This, of course, 
is what has been observed. 


Assuming that the pressure returns to the value just in front of the 
shock, а simple calculation shows that the turn of the flow through а 
Prandtl-Meyer expansion has to be 


-в( (o) š FIO) 


where 1002) (0) and sult) (0) are, respectively, the velocity ull) 


just before and after the shock. By means of solutions (27) and (51), 
the turn required is е. Now the direction of the flow before the expan- 
sion is є(1 - BA); therefore the total inclination of the velocity 
vector at OM is 


е(2 - ВА) 


A solution in region 3 subject to these conditions is found to be 
(&ppendix B): 


(1) A oA (x-By) 


Po = Ge + By) - (кх - Ву) + X (52) 


% 
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The velocity in region 3 cam then be given by 


sult) = Ae» (х-Ву) 
тї) a Ys. 8Ae^ (x-By ) 


Were the perturbation velocity to remain finite at infinity, X would 
have to be negative in region 3, as x - By 2 0. 


Viscous Solutions 


As potential solution (52) in region 3 is a superposition of two 
different types of waves, namely, step waves and transmitted waves, it 
is expected that the viscous solution will also be composed of two dif- 
ferent parts, each associated with a special wave in the external field. 
Since the origin and character of the disturbances are different, they 
can best be discussed separately as the following: 


Im view of the fact that the step waves are opposite in sign, they 
do not change the pressure but give rise to a uniform vertical velocity 
in the potential field. The flow in the viscous layer due to this uni- 
form deflection may vary rapidly in the y-direction but certainly not 
in the x-direction because of the constancy of pressure. The main 
effect of the step waves, first of all, then, is to produce a constant 


deflection along the edge of the viscous layer. If the velocities due 
to step waves are u UJ and v, V3, the problem should be solved 


subject to the boundary conditions: 


ug (1) = Жез = 0 when y = 0 
u, (1) = 0, Жез = 2 at outer edge of viscous layer (53) 
u, V) = 0 when x = 0, y£20 


It must be borne in mind that generally the flow in region 1 has 
been separated, with a considerable region of backflow in the neighbor- 
hood of the point of incidence; the resultant flow in the incompressible 


layer must be small. If u 1) and. y, V are considered as additional 


— n 0 M HÀ € 
2 ас p c we ee 
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perturbations due to the step waves, the problem сап then be simplified 
by neglecting the inertia forces in the,neighborhood of the plate. Тһе 


perturbation stream function y, (1) then satisfies 
A ay VH -0 (54) 
in the incompressible layer OS у & Yi. 


Ав the objective here is to demonstrate approximately the effect 
of the step waves on the viscous layer, only simple solutions will be 
considered. А simple solution that satisfies both the first and third 
of conditions (53) is clearly 


for y in O Sy S y}, where D is а constant. The velocities are 
thereby 


(56) 


Accordingly, the pressure, or the pressure gradients, in the incompress- 


ible layer due to the step waves is of the order К- and for the pres- 
ent approximation will be neglected. I 


In the compressible layer where the viscous effect ів not as impor- 
tant as the compressibility effect, a flow that does not associate with 
& pressure rise can either be & uniform field or the one with & vertical 
velocity varying linearly as y. According to equation (37), the fol- 
lowing,' with the second of equations (53) satisfied, is & solution for 
A < y < Yor Yo being defined as the outer edge of the compressible 


layer: 
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(57) 
vg) 


The fact that equation (57) is a solution shows the importance of both 
the compressibility and the viscosity in this layer, which allow the 


independent variation of oC? and т(1), The joining condition at 


the interface gives the constant D a value -6y3 /Yo: 


Thus, it is seen that the deflection of the flow due to step-waves, 
though it contributes no pressure, induces a forward velocity in the 
incompressible layer. Since it increases linearly with x, the back- 
flow should be expected to be reduced in the downstream direction. Of 
course, the possibility exists that sn additional pressure might also 
enter if the inertia forces, though small, were not neglected.  Never- 
theless, it can be safely stated that the pressure is always of the 
secondary importance in this case and hence the above conclusions will 
remain valid. 


On the other hand, the perturbation velocities u 1) and y, (4) 


due to the transmitted wave, according to equation (52), are subject to 
boundary conditions 


u, VJ = v (1) 0 when y = 0 
a V) Ae^* v, V) = -pAe^* when у = х (58) 
(1) _ ,(3) when х= 0, у20 


1 


To solve this problem, the viscous layer is again subdivided into incom- 
pressible and compressible layers. Since the boundary conditions and 


the differential equations are the same, a VJ and y, VJ will have 


the same forms ав those given in equations (33) and (hh), with eigen- 
value satisfying equation (48). Now if the perturbation is required 
to vanish at positive infinity, X should be negative. But were A 
negative, solutions (33) and (lh) would be highly oscillatory, as the 


ы 


У 
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arguments of both Bessel and the exponential functions involve & fac- 
tor VAR. That is, the solution of the incompressible layer would 


| . 3/2 
involve s (2p, oar y J and that of the compressible layer would 
3 


have cos VAR y and sin VAR y. For large Reynolds number, these types 
of solutions for a steady laminar flow must be rejected as impossible. 
By the same reason, complex roots are excluded. The only alternative is 
to accept the positive А. 


If А is positive, it follows that the velocity guy and pres- 


gure p, 4) will continue to increase until the process of linearization 


breaks down. This would imply that in the potential field the train of 
compression waves, immediately following the expansion, will grow expo- 
nentially with x. The flow direction as well as the curvature of the 
streamlines will also increase sharply. For such a flow, it is well- 
known that the Mach waves will converge and form an envelope.  There- 
fore, in region 3 & shock eventually must be developed. On the other 


hand, in the viscous field, as y V) + u, Q3 for any constant y will 


increase with x, backflow, though slightly reduced by the step waves, 
will become stronger in the downstream direction. Now, it has been well- 
established that the laminar velocity profile with & point of inflection 
is highly unstable. As the pressure continues to grow at a Reynolds num- 
ber generally above the critical value, transition must occur after a 
critical pressure gradient is reached. The flow from there on cannot 

be theoretically studied without considering unsteady flow. 


It is therefore concluded that in the case of an incident compres- 
sion shock, laminar flow is not possible for the whole viscous layer 
and transition always occurs. This appears to be in complete agreement 
with the present available experimental observations. 


Character of Flow After Transition 


Although the flow from a certain point on is unknown, the inter- 
esting fact is, however, that the flow up to the point of transition is 
very insensitive to what happens beyond the point of transition. Іп 
the previous sections, the problem has been reduced to depend only on 
one constant A on which the quantitative behavior, but not the char- 
acter of the flow, depends. А quite similar conclusion has also been 
reached in reference 1 and seems to be an experimental fact (refer- 
ence 7). In one of his experiments, Liepmann introduced an expansion 
wave immediately after the incident shock, and the observed upstream 
flow field was practically unchanged. Therefore, in order to account 
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for the observed pressure distribution, say it is sufficient Just to 
take approximately the effect of the flow downstream of the point of 
transition. 


Assume that the second shock starts at the edge of the viscous 
layer апа, by linear theory, has the same strength в.  Inasmuch as the 
direction of the flow on OM is constant and equal to «(2 - BA), if 
the flow after the shock TR (fig. 2) is parallel to the plate, the ` 
direction in front of the shock at large distance must be є. This con- 


dition gives A= B Moreover, as the velocity vector is turning away 
from the plate, the pressure continues to rise and will be stopped only 
by transition. Іп that event, the sudden increase of shear of the vis- 
cous layer will thicken the viscous layer. This thickening will make 
the flow expand and consequently a drop of pressure will ensue. After 
this, the streamlines will gradually level off and the reattachment of 
the partially turbulent layer, if not yet &ccomplished prior to the 
transition point, will be certain to follow. If the point of transition 
41s chosen to coincide with the point where v = O along the boundary 
streamline, the location is then determined by Ax, = 106,2. 


For values of x greater than x), the flow in the immediate 


neighborhood of transition would have a greater influence than that far 
downstream where the flow is more uniform. Since the flow in the tran- 
sition region is of a boundary-layer type, namely, the backflow ceases 
to be a factor, it can again be approximately represented by the inte- 
grated effects, such as the momentum and pressure. If the flows before 
and after the point of transition have the same pressure and total 
momentum, the dynamical equilibrium can then be maintained. According 
to the solution given in reference 1, the pressure eon in the 
transition region is approximately 


Ах x 
іе” + Ke? 
where 24 and А are negative constants, being functions of M 
and R and Ку and Ко are integration constants. By the conditions 
that at x = ху the pressure and its derivative are continuous, Ki 


and Ko can then be determined. 


It must be emphasized again that the conditions stated in this sec- 
tion are tentative. No obvious reasons beyond the ones outlined at the 
beginning of the section can be given at this moment. 


——— Á --- -— ---- ——————  — — —— HY. = — = 
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NUMERICAL EXAMPLES 


Ав a numerical example, the case of М = 2 and R = T[h is pre- 
sented. First, the determinantal equation (50) is solved numerically 
for only one Mach number Z but at different Reynolds numbers. The 
results are shown in figure 3. At this Reynolds number, À is 0.0467. 
The velocity u(x,y) апа the surface pressure р(х,0) for this X are 
calculated. It is seen that for the shock strength of є = -19, the 
boundary layer first separates at Хх = 1.106 (fig. 4(a)) and, subse- 


quently, backflow sets in. When the shock strength is increased to 
€ = -39, say, the separation occurs at a much earlier station, namely 


at Xo) = -22.4. Тһе region of backflow is proportionally wider 


(fig. 4(b)). However, if the step waves are taken into account, then 
in the case of є = -1° there would be no separation (dashed curve in 
fig. Қа) for уу = 3 and y, = 9); whereas in the case of 6 = -3? 


the flow at the same location remains separated but the backflow has 
been very much reduced (dashed curve in fig. h(b)). Теве results show 
quite the same characteristics as the experimental measurements by 
Ackeret, Feldman, and Rott (reference 8) of the velocity distribution 
over a curved plate. 


The pressure distributions over the plate for these cases are shown 
in figure 5. For the weaker shock 6 = -1°, surprisingly good agreement 
with experiment (reference 6) is obtained. This very fact seems to 
justify the present assumptions regarding the structure of the viscous 
layer. In the case of stronger shock, for example, є = -3°, there is, 
however, a distinct difference between theory and experiment. Theoreti- 
cally, the pressure would still decay exponentially upstream but would 
begin with a larger amplitude. Experimentally, it was found, strangely 
enough, that, over a considerable range of the upstream disturbed region, 
the pressure first decreases very slowly and then decays more or less 
exponentially. This "pressure bump" seems to be characteristic of the 
pressure distribution in the case of interaction between stronger shocks 
and the laminar boundary layer. This bumpy character in the pressure 
distribution, by all evidences, must be attributed to the nonlinear 
effect of the flow. This will be exhibited in the following section. 


APPRAISAL OF HIGHER-ORDER EFFECTS 


w 


It has been shown that, if there is no backflow in region 1, the 
pressure disturbances will decay exponentially. When the shock strength 
increases, however, the pressure in the disturbed region becomes much 
higher and drops much more slowly than that predicted by the theory. 
This appears to be due to the fact that when backflow develops, there 


“е m — — - — — F. н н — = — —— — .— r — sn — 
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will be an underestimate of the perturbed velocity and, consequently, 
a much lower pressure. 


To estimate the effect due to higher-order terms, а second-order 
solution is found and is given in appendix C. According to the second- 
order solution, the correction terms for wall pressure would be 


Руел + Baer 


In view of the fact that there is an overestimation of the first-order 


pressure by the condition А = B =, the correction term is expected to 
be negative for small values of x апа positive far upstream. Then 


Pa > O and P, < O. Since a AE increases faster than x^ and eX, 


when x 15 positive, the correction would be negatively large far down- 
stream. It has been pointed out previously that the iteration process 
fails for positive values of x; this solution shows further that it 
diverges oscillatorily. Tn order that the solution will behave properly, 
it will have to be carried to an odd order, such as one or three. To 
third order, the correction for wall pressure would be 


P ехх + p; ex 3Ax 


1 + ы 


If Ру = 0.15, Р» = -1.83, and P3 = 1, the wall-pressure distribution 


for € = -39 would resemble the curve as shown in figure 6, which does 
exhibit the same character as measured by experiments. This, of course, 
cannot be considered as conclusive evidence, but at least it shows that 
higher-order terms such as given &bove do have the possibility of 
&ccounting for the observed bebavior of the wall-pressure distribution. 
It is very important tbat this step be carried out. 


DISCUSSION 


This study being intended as sn exploratory study, the numerical 
results obtained are not expected to be exact, but only accurate enough 
to insure the conclusions. Ав far as the first-order solution is con- 
cerned, the entire problem depends on the determination of one param- 
eter А, which was calculated on the basis of two main approximations: 


co 
(a) the contribution of the integral / 1922 (9/2 аз is negligible 
74 3 
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and (b) hy, ів small such that terms of order (Xy, y^ and higher can 


be neglected. Тһе subsequent solution for X satisfies this criterion. 
The question remains, however, whether Лу; = n(1) is also a solution 


" 
or not. If Ay, = n(1), the integrals J + 2/2 at Ky (Š 13/2) an 
0 3 


would require а more elaborate evaluation. Тһе possibility exists ала 
is worth testing. 


As from the outset the boundary-layer approximation is rejected, 
it is of interest to examine the character of the perturbed flow in the 


light of boundary-layer theory. According to the boundary-layer &pproxi- 
-1 
mation, if the velocity ш! = o(U, ) and Ә/дх' ш olz )n then the 


velocity v' = ouR!) and o/ày' = o(r-1g). Now, from the character 
of the solution for perturbed flow in the present theory, д/ àx ! = о(т-1н3/%) 


and Ә/ду' = о(тг1ң°/%). By comparison, the derivatives for the perturbed 
flow are much larger than the corresponding derivatives of boundary-layer 
theory. Moreover, from the equation of continuity the ratio GA G) 


of the perturbed velocities is о\в-1/2). Ав in the vicinity of the 
wall 4) a 60) y(t) is mich larger than the vertical component in 
the boundary layer. Since in this case 0/dx' < д/ду! and y GJ) <«(1), 
a set of equations analogous to those for the boundsry-layer flow can be 
deduced as long as backflow does not take place. Therefore, for expan- 
sion and even weak compression incident waves, a much simpler problem 
would be feasible. 


Finally, it might also be noted that, owing to the assumption that 
the flow far away from the plate is inviscid and irrotational, there is 
introduced a sharp discontinuity in higher derivatives at the demarcation 
line between the two flow fields. Because of the presence of shock, how- 
ever, a first-order discontinuity is also expected, because, by the 
assumption of frictionless flow at large distance, a smooth transition 
from small shock thickness at large distance to a larger one at the 
vicinity of the wall is precluded. The viscous-layer concept then 
idealizes the situation by taking the shock as a discontinuity in the 
outer field but continuous in the viscous flow. The effect of this is 
exhibited in a discontinuity in slope of the streamlines. This picture 
is entirely in agreement with the observed flow patterns. 


According to the numerical example, the step waves tend to weaken 
the backflow downstream of the point of incidence but are unable to make 


ме - — --- — --- 2 —— € =) —— — = — — 
- ---- nm  ——  — — A LR E Tate RS UT ER m nre I ag  — —À ——— — OÚ a — ——— тейін 
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the flow reattach if the separation has occurred upstream of the point 
of incidence. This may very well be due to the poor approximation 
inherent in the Stokes flow which шау underestimate the rate of change 
in the y-direction. If & more exact solution is given, reattachment of 
the separated flow might be accomplished in the laminar regime. This 
point should be considered as open. 


CONCLUSIONS 


The following conclusions are drawn from an analysis by the 
differential-equation method of the interaction of an oblique shock 
wave with a laminar boundary layer along a flat plate: l 


l. The pressure perturbation decays exponentially forward of the 
point of incidence апа the distance of pressure propagation varies with 


2 


different Масһ and Reynolds numbers as p 3/9 -3/5 ( where P is Mê - 1 
M is Mach number, and R is Reynolds number ) and very slowly with the 
shock strength. 


4. If the shock strength is strong enough, separation of the flow 
always occurs. For given Mach and Reynolds numbers, the separation point 
depends strongly on the shock strength. 


3s The pressure &t the point of separation varies with Mach and 
Reynolds numbers as g^ 1/29-1/2. 


h. In the viscous layer, laminar flow is not possible everywhere, 
no matter whether the incident shock is strong or weak. Іп the distance 
of about two or three boundary-layer thicknesses, transition would occur. 


S. The curvature of the streamlines after the shock is positive and 
the Mach waves in the potential field must coalesce to form a shock which 
approaches asymptotically the regularly reflected shock in the inviscid 
fluid. As its position depends on the point of transition, the exact 
location cannot be predicted by the present theory. 


6. From the calculated pressure distribution over the wall, it is 
definitely proved that the observed overcompression of the wall pres- 
sure is a consequence of the positive curvature of the streamlines and 
the expansion is associated with transition. 


T. The observed "bump" in the pressure distribution in region 1 for 
strong shocks is definitely a nonlinear effect and is am indirect conse- 
quence of separation. | 


Cornell University 
Ithaca, N. Y., January 11, 1952 


GE 
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APPENDIX А 
PERTURBATION ОЕ AN OBLIQUE SHOCK 


let Up yp Pp and 04 denote, respectively, the parallel 


and normal velocities with respect to ап oblique shock, the pressure 
and density in front of the shock, and, by replacing 1 by 2, the corre- 
sponding quantities behind the shock. These variables are related to 
each other by the Rankine-Hugoniot conditions: 


P1 19, = P2 oy 


(a1) 


If the flow in front of shock is given, equations (Al) yield the fol- 
lowing solutions: 


y-lt* 43-) (a2) 


TOM z - —— E — — пн e m 
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where aj is the speed of sound.in front of the shock. If a, and U 
are the speed of sound and velocity at infinity, then а can be 
expressed in terms of velocities by the relation: 


со 


2 - if 2 2 2 = г 
aj — (us tv?) = P +2 Uo 
Furthermore, if œ is the angle between the velocity (mv) amd the 


shock, it is easy to show that the Cartesian cómponents before and after 
the shock are related by 


2 
" 2 2 1 
uo = 1 + — (ua sin o + Үз sin o cos o| ——— - 1 
44 gin o 
(АЗ) 
2 a,” 
Vo = 1 - (а) сов о sin o - V4 sinfo)| 2 - 1 
1 2 2 
44 gin 0 


with 12 = uj + 7,2. 


In the case of weak shock, the shock angle œ сап be expanded in 
terms of une deflection angle є. Then cos o and віп 0 are shown 
to be 


Eee s. 


(а) 


(о +1) al + G2 му. „к. к 
hg? 8g? 
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Now since the velocity components uj and уу can be expanded: 


1) T cu, (2) + 


= 0) a OD 4 


Uy = ] + cu, Í 


Уі 
If o has the expansion: 


CONTENT OUR 


and 209) is defined by equations (Ан), then a кес reduction 
gives 
uo = l + cu, 19 + eu?) + 
оле eg + VOLL. (45) 
where 
(1). h _ X (7-1 (1) ‚1 1p2q(1) 
а : Du às deem 
vo (1) „y (2) _ ue "E - 1 5) у 38292) 
°: 1 y + 1 l 2 E d Tos d 
(1) 
(2).]3. + 0-1,1, (2) 2 isa a) s ) 
™ з = a EJE dE m B 
Jg i» lh 1), (1 
S zer QQY. (v, y PB ss (1), 
(46) 
b 1(7+1м№ _ 1) - (м2 - 2) ot) 4 ET 2 = 20 - 2)(,, (1))° 
p2\ og 82 (y 1M (у +1) 


v (2) =- Se 
li 2 


арғы (2) 
(y + (y + 1)⁄2 


ur zuo 10) « SEE a) a 


Say @) + vy (2) + 28 eer + 2M (a)? - (01 


PRORA ‚у-у G) уз, wea? | ca) _ 
В (y + 1)м 2 В (y + 1)№| > 


, 2В(1 + ки 
(y + 1)м 
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APPENDIX B 
POTENTIAL FLOW IN REGION 2 


If the flow in region 2 is irrotational, there exists & poten- 
tial 29(х,у) and, if the shock conditions are satisfied, it must 


possess the expansion: 


oP =X + € ,90) £c ,9 0) Че op. 4x (B1) 


For the first-order solution of the shock conditions, namely, on & line 
x + Ву = 0, 


,9, 0 = nue sue A 
— (у + 1) 


(B2) 


(1). , (1) _ 4B е -1. DE ЕТ қаға) 
(у + 1)м 


Now on the shock 
v, Q) = -BAe?^x 


If the perturbed shock angle is 


wit) = geMx-By) 


the solution 292) must be 


. 2 


gt) = Ha + Ву) + te^ Py) | (вз) 


+ 
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Substituting nO and 2% іп equation (B2), the two equations 
yield uniquely | 


(Bh) 


The solution in region 2 is thus & superposition of the transmitted 


waves od (x-By) and a uniform step wave x + By introduced by external 
agency. 


Since this solution will terminate on the Mach line OM, the solu- 
tion in region 3 will be found by the condition that on x - By = 0, 
(172.42) 
qu у (В5) 
By this condition, the Prandtl-Meyer flow requires а deflection 


59 = -ep yu (0) - „ч (о) (B6) 


Solution (АЗ) gives 59 = €. Consequently, the initial direction of the 
flow in region 3, to the first order, is 


„у (0) + 5% = (2 - ВА)є (B7) 


1-- -- —-. i Na nF — -— 
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APPENDIX C 


SECOND-ORDER SOLUTION 


Potential Field 


If the first-order potential is of the form: 


pU А eà (x-By ) 


it can be shown that the second-order equation (19) admits a solution 
which reduces, at the edge of viscous layer y = O, to 


ul?) = Beò” + ge Ax 


(c1) 


2 
vie a. И ПР ы Әм QA 
hp 


2 
(y + MA „2\х 
hg 
2 
(у + LM AS ny (x-By), 
hp 
2\х Ах 
reason that both е and e are included in the complementary part 


of the solution is that the solution in region 2, according to equa- 
tion (AG), contains both solutions. 


where B and C are arbitrary constants and the term - 


is derived from the particular invegral - The 


Incompressible Layer 


Assuming that, to the second order, the density in the inner layer 
remains constant, the system of equations (15) simplifies to 


A луб) - бу Ax, 2) = TAR (1) - AA) (C2) 


x 
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where the second-order stream function y 2) is defined as 

42). y @) 

ge) 2 y (2) 


By substituting y(1) into the right-hand side of equation (C2), the 
particular integral can be found; this is 


ғ... 7 ды 
a = = | (cos тп)2(%) ат e 
* (28. )#/3 0 


ww 
сов Th 


ў 
ИК 6 (sin 79)2 (9) аў |е? (C3) 
тал 2/3 D 


where the function 7(П) stands for 


(4) = 2 4/22 99/2) | i qx (2 93/2) a(n) asi] - 
3 3 


ағала ты tn 


—— ——— ——— -- — ——— F wn — 07 


———  —a,sF— D — —— — — — r — w == — — — — = 
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with 


2 
A(n) RS Lean)? - п/к (2 13/2) = m) 51/2 (cos TT) ку ат + 
3 3 


Бо с 


(sin "0 10/2 (sin тт) Жү 1 - 
3 


7 
i (2 13/2) Z Tn) J, „1/2 (cos тт Жү dy - 
3 3 
(сов mf 11/2 (sin L5 (C5) 
0 3 


and the variables: 
7 = (ғал) З(у ғ 24 | (сб 
) 


Y = ру, Aa) (cT) 


The general solution, consequently, сап be written as 


92V 0, v, (c8) 


where the complementary integral "E ) ; by the first-order solution, 
will be given by 


Vo? - сор) (пех + c, y (е (c9) 


----- -- —  . — 
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I 
Since the form of VU) is known, and, furthermore, q (o) = "n (0) = O. 
the solution y (2) is determined with two arbitrary constants C3 and Су. 


With y (2) determined, the second-order pressure 5) (х,о) on the 
plate can be written, similarly, &s 
р(2)(х,о) = p eà + рге ХХ t рге Хх (C10) 


where 


__ (ға. ) / 3342/3 г E 


pois (2) Gr) Jo EE 3/2 JAn) а 


and P4 and Pj take the same form as the symbol P for the first-order 


pressure p (1). 


Compressible Layer 


In the compressible layer, p 9? = (0) = m(0) = l. Again, if 
2) 


н! = 0, equations (15) become 


(Ú 2) . 2), -24(,9 (20) 090090, «С 309), 


vagy (1) Ж ullar (1) 


"MOORE OE: ny (2)) - -(1 + aÊ) G) + 


Еж (1-2 АРУЫ QA y(t) - ar(r +2 g a = laf a) 


му - x40) 1G, QD . M" 


(c11) 
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2) 2) 


Similarly, by introducing p and с through the relations 


| u (2) = ғ, (2) + a, (2) 


(2) в (0) a @) 


y 


Gy 


there result the following equations: 


y (Pag, 2) Dag, (2) 


AOG, " 6,02) sAn) a зе |6) " a, (2)) " Pag В 


(C12) 
(1+ эё), @ Q2 , E ит ле), G) _ 
(+ ZE Ba a G - ape? а) аа), Q) 
The solutions p) and (2) can be written in an analogous 
manner. "That is: 
ШЕ)» (2) p (2) 
(C13) 
2). 5002) + a (2) 


Here JAG, and со 2) аге known from the first-order solution and 
ғы 2) ала ср 2) are given ав follows: 


F (2) = Can Koe 97 + ке 2Y) АЖ 


(C14) 
-су+2Ах 


= Kye 


ee  — - — 
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ol 


where the constants K аге defined by 


2+ ( 1+ (82/62). 2) |+ (y - axes? + (у - 1) + 2 +20," 2, 
“а(ь2- 62) Ba”) 8 y, 


l+3 -0 


1 - (8/82) Ë - (7- 1) |8,2 +2 + (y - UM a 
By i S 


1 - (р/н) Е №2 + (AB, - 7 
K, = ——— [x84 - с} |2@ - у 
st to eS 


- anma — — —  - ee i 


(22,2 - DE - xl + polk + eal 


alı - zl) + (xP, - sy t 2(1 - юа?) - (хв) (8, - e) 


+ 


g&1 + > - = ioco 19 + Ы o - ane es 


(= - 862 T 20 + 2) + adı + ум + Male - 1) - 4) x 


CoA 
2 


—TV + 
(AB, +a) - n Nd 


1 + а/ы) Y. + 24-2 + олуы? + (Ву - е - AS — - š - Ex + ға) | 
) 


2 + 2 ж. 82(1 - 2) + 2(X84 4 “Қ - =) - (XR) 138, 4 e) 


SG + në) - SE rp o - ay - A] +a + zty -aël + 


SACRE DRE E 


CoA 


(X81 + сў - ap 


— — -- — ——s a — — — -- жу ле ренин 
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-1 _ 22р? 
на) (y - 1) ей J- 0-2 p 


Вус 


G5 - о®)(1 + в/в) 5 uae (x ú l 


i terri te шк Lc p Mu e E EEN E E E l u тк 
у anfaz - 28-1) + 208, + cJ. - wart) - (x)7168, + 2d 


(22,2 - s^) - в/в.) к а Se + 2) 


ӘКЕ - art) + 2(XB. - sy 2 вті) ов) оз, _ sy 


and the constants K and В. stand for 


8 ум (C15) 


The second-order solution now involves six arbitrary constants and, by 
the conditions at the interface, all but one can be determined. 


* 
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Figure l.- Diagram of flow field. 
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Figure 2.- Sketch of flow field according to first-order theory. . 
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Figure 3.- Curve of АХ against R for М > 2. 
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(а) в -19; M = 2; R = 774. 


^. Figure 4.- Velocity distribution. 
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(b) є=-3°;М=2; R = TTM. 


Figure 4.- Concluded. 
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Figure 5.- Pressure distribution on plate for 6 = -3° and є = -19, 
M = 2; R = 774. Experimental values were taken from reference 6. 
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Figure 6.- Pressure distribution on ушы for third- апа first-order 
terms for є = -39. M = 2; R = 774. 
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